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Discovering baryon-number violating neutralino decays at the LHC 
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Recently there has been much interest in the use of single-jet mass and jet substructure to identify 
boosted particles decaying hadronically at the LHC. We develop these ideas to address the chal- 
lenging case of a neutralino decaying to three quarks in models with baryonic violation of R-parity. 
These decays have previously been found to be swamped by QCD backgrounds. We demonstrate 
for the first time that such a decay might be observed directly at the LHC with high significance, 
by exploiting characteristics of the scales at which its composite jet breaks up into subjets. 



The LHC's scope for the discovery of supersymmetry 
broken at the TeV scale [IHE] has generated much inter- 
est. Certainly, the potential prize is great: TeV-scale su- 
persymmetry could solve several puzzling problems and 
answer a number of open questions in modern particle 
physics, such as the fine-tuning of the Higgs mass, the 
unification of forces at high energies and the nature of 
dark matter. Effort has mainly been concentrated on 
investigations into the discovery reach and possibility 
of parameter measurements in the Minimal Supersym- 
metric Standard Model (MSSM) and various more con- 
strained versions featuring a weakly-interacting and sta- 
ble Lightest Supersymmetric Particle (LSP), which gives 
a missing-energy signature. Candidates for the LSP in- 
clude the lightest neutralino, Xii and the gravitino. 

However, the gauge symmetries of the Standard Model 
(SM) also allow for dimension-four terms in the superpo- 
tential of the forms 

XijkLiLjEk + \ij k LiQjDk + \"j k UiDjDk, 

which violate the R-parity that is imposed in the MSSM. 
Non-zero values for the couplings A could imply drasti- 
cally different phenomenologies for supersymmetry at the 
LHC, allowing for the decay of the LSP, or any heavier 
sparticle, directly to SM particles. 

Broadly speaking, we can classify the RPV models in 
terms of the dominant coupling A and the identity of 
the Next-to-Lightest Supersymmetric Particle (NLSP), 
assuming that a gravitino is the LSP. The lepton-number- 
violating (LV) couplings A and A' generally give MSSM- 
like signatures with missing energy from neutrinos and/or 
extra lcptons in the decays of the NLSP. The exception 
is a slepton/sneutrino NLSP that decays into two quarks 
via a dominant A' coupling. Scenarios such as these 
should be easy to extract from the SM backgrounds, as 
shown for a neutralino NLSP in Q. 

However, dominant baryon-number-violating (BV) 
couplings A" are more difficult to deal with, due to 
the large hadronic activity expected at the LHC, which 
threatens to drown decays such as Xi — > QQQ- The QCD 
background for jets with px > 500 GeV and a Xi mass of 
0(100 GeV) is about two orders of magnitude higher 



than the signal. The background's non-trivial shape 
means that it is hard to establish whether a small devi- 
ation from the expected background is a signal of some- 
thing new, or simply a defect in one's understanding of 
the background. Some success has been reported by re- 
lying on the production of a high-p-r lepton in the decay 
chain leading to the NLSP [6-8], but ideally one would 
wish to demonstrate the feasibility of signal isolation and 
mass measurement in a less model-dependent manner. 
Otherwise, the fear is that supersymmetry could escape 
discovery at the LHC by cloaking itself in BV decays. 

In this Letter we investigate that problem by looking 
for jets from the decays of very boosted sparticles via BV 
couplings. Such decays give rise to composite jets made 
up of two or more collimated subjets, with a jet mass re- 
lated to that of the original sparticle, with specific prop- 
erties predicted for the scale at which the main jet sep- 
arates into subjets. Similar techniques have previously 
been used by the authors for analysing WW scatter- 
ing [9l ITol ) fo r detecting massive boson decays in MSSM 
scenarios [111 and in Higgs searches In addition, 

a number of other techniques for separating hadronic 
decays of heavy particles from QCD backgrounds have 
been suggested by other groups [13H19I ]. While some of 
these studies are applicable to general searches, they all 
worked with examples of particles with known or well- 
constrained mass. The present article tackles a new 
problem: how to identify a hadronic resonance of un- 
known mass in an explicitly scale-invariant manner. The 
techniques presented here in a supersymmetric scenario 
with RPV clearly have applications to any hadronically- 
decaying massive-particle resonance that can be pro- 
duced far above threshold, and are promising for broad 
use in the challenging LHC searches for hadronic decays 
of new particles. 

We focus our investigation on the CMSSM benchmark 
point SPSla. We also look at six other CMSSM points 
with larger sparticle masses and lower cross sections, 
that lie along the corresponding benchmark line [20| . 
These encompass sparticle masses from (m^o ,niq,mg) ~ 
(96,530,600) GeV to (161,815,915) GeV, and cross sec- 
tions from 47 pb to 4.6 pb at leading order. We note that 
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the gluino and squark masses are much larger than the 
neutralino mass, and can hence yield a highly-boosted 
neutralino in their decays. 

R-parity violation is incorporated by setting the domi- 
nant coupling to be A" 12 = 0.001. This coupling is chosen 
to be difficult: no heavy flavours are present to help tag 
the correct jets and the coupling is chosen to be relati vely 
large, so that decays do not lead to displaced vertices [21j . 

In order to simulate sparticle pair-production events 
at the LHC with RPV decays, we use the Herwig 6.510 
Monte Carlo event generator [^IJp] with CTEQ 6L 
PDFs, and use Jimmy 4.31 [27t| for the simulation of 
multiple interactions [28[ . This is interfaced to the Fast- 
Jet 2.4.0 [HIU jet-finder package using the Rivet [H| 
framework. Our background sample, consisting of QCD 
n-jet events, tt, W+jet, Z+jet and WW/WZ/ZZ pro- 
duction is simulated with the same setup, with the excep- 
tion of the n-jet events where we use Alpgen 2.13 (32J 
interfaced to Herwig / Jimmy to generate up to 4-jet fi- 
nal states with the MLM jet-parton matching scheme. 
For the QCD event generation we also require ~^2pt > 
600 GeV to reduce the large sample size, where the sum 
is over all jets. The leading-logarithmic parton shower 
approximation that is used has been shown to model jet 
substructure well in a wide variety of processes [33j-|38j. 

For both signal and background we generate a number 
of events equivalent to 1 fb _1 of LHC data at 14 TeV 
CM energy. No attempt is made at detector simulation 
through finite calorimeter granularity, but we do impose 
a geometrical acceptance cut on jets of 1 77 1 < 2.5. 

To illustrate the variety of approaches that can be 
taken with subjet studies, we consider two complemen- 
tary analyses: the first is an exclusive analysis with 
the algorithm 39, H(| and searches for substruc- 
ture in two jets, one for each neutralino expected in an 
event. The other is an inclusive analysis, with the Cam- 
bridge/Aachen (C/A) algorithm [Il|,[42j], which examines 
the substructure of just the hardest jet. We find that the 
corresponding inclusive (exclusive) analysis with the kx 
(C/A) algorithm give similar results. 

The kx algorithm defines distances dki = 
min (p Tkl p Tl )(ARlJ R 2 ), d k B = P Tk , and sequen- 
tially merges the pair of objects k, I with smallest du, 
unless there is a smaller dkB, in which case k becomes 
a jet. The constant R sets the angular reach of the 
jets. Since the kx distance is just the relative transverse 
momentum between objects, the mergers of interest for 
a decayed heavy-particle tend to be the last ones. This 
was exploited in [9T4llj|. where a dimensionful cut was 
placed on the dki scale of the last merging in the jet, d\, 
in order to preferentially select boosted W bosons over 
QCD jets, which, for a given mass, have smaller d\. 

However, our case differs from 0-[ll| in two respects. 
First, we are searching for an object of unknown mass, 
which means that we should avoid biasing the search with 
a dimensionful substructure cut. A good alternative is to 
cut on a dimensionless variable normalised to the jet mass 
m 3 , yi = diR 2 /m 2 . Secondly, the neutralino has a three- 




FIG. 1: The y-values from the &t algorithm with R — 0.4 for 
the last and next-to-last merging for QCD jets (left), and jets 
matched to neutralinos (right). Also shown is the proposed 
cut line. Both distributions are normalized to unity. 



body decay, in contrast to a W-boson's two-body decay. 
This suggests that one should cut on the properties of the 
last two mergers, i.e. on yi and 2/2 ■ Their distributions 
are illustrated in Fig. [TJ The left panel shows QCD jets 
with px > 300 GeV, while the right panel shows jets 
matched to neutralinos for the SPSla benchmark point, 
the jet being within AR < 0.3 of the neutralino with 
p T > 300 GeV and mass 90 GeV < m,j < 120 GeV. Even 
after the hard cut on jet pt, there are clear differences 
between neutralino jets and the QCD background. 

For our analysis with the kx algorithm, we take R = 
0.4 and use the following two cuts: i) at least four jets 
with p T > 300, 300, 100, 100 GeV and two of the high-p r 
jets (px > 300 GeV) should have significant substruc- 
ture, y2 > — 0.13yi + 0.05, consistent with fully or par- 
tially contained neutralinos. The choice of R is a compro- 
mise between capturing sufficient signal, favouring large 
R, and not smearing the mass peak with particles from 
the underlying event, favouring small R [431 ] . Cuts on a 
third and fourth jet are motivated by the expected pres- 
ence of two jets from the squark or gluino decay. We 
have verified that our simulation of QCD events with ad- 
ditional jets is consistent with NLO calculations in [44j . 
The requirement of high px should ensure high trigger 
efficiency, as well as good collimation of the neutralino 
jets. 

The mass distribution for the neutralino candidate jets 
after cuts is shown in Fig. [2] (left). The x? is clearly vis- 
ible as a perturbation on the rapidly falling QCD back- 
ground. This analysis reconstructs 3.3% of all neutralinos 
with pt > 300 GeV (constituting 17.8% of neutralinos 
produced) in a 20 GeV mass window around the nominal 
mass; most reconstructed neutralinos have a pt some- 
what above the 300 GeV cut. For the QCD jets the 
analysis accepts 0.023% of all jets with pt > 300 GeV in 
the same window. 

Even though we have introduced only dimensionless 
jet substructure cuts in the above analysis, the back- 
ground distribution also has a peak near the neutralino 
mass. This is a consequence of higher-order perturbative 
effects [13, |45| and the peak position is determined by 
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FIG. 2: Jet mass distribution (points with error bars) for the 
kr algorithm (left) and C/A algorithm (right) after all cuts. 
Also shown are the contributions from QCD (black), super- 
symmetric events (red) and other SM backgrounds (green). 




FIG. 3: Estimated sensitivity for 1 fb _1 as a function of \i 
mass for various choices of jet algorithm and size R (left), jet 
mass distribution for a squark search using the C /A algorithm 
with R = 0.5 and background estimation by sidebands (right). 



their interplay with R, the jet pt and substructure cuts. 
Increasing R, which improves the signal reconstruction 
efficiency, causes the background peak to move to larger 
mass, potentially directly beneath the signal peak 

To avoid this issue, we consider the C/A algorithm, 
which successively recombines the pair of objects closest 
in ARki, until all objects are separated by more than R, 
at which point they are the jets. Because the ordering of 
C/A mergers knows nothing about the momentum scales 
involved, one cannot rely on the properties of the last 
mergers to tag relevant substructure. 

Instead we recurse through all mergers storing those 
subjets that are sufficiently symmetric, having z = 
mm{p Tk ,p T i)/(p Tk +pti) > Zmin, ignoring in the recur- 
sion the softer of the two subjets when z < z m in- From 
this we identify the two mergers that have the largest 
jADE-type distance, d kl = pTkPTiAR\ Xl which is related 
to mL if k and I are massless. If the subjet resulting from 
the merger with smaller d J (labeling this subjet "be") is 
contained within the subjet from the other merger (la- 
beled "abc"), i.e. the subjet "abc" is the result of further 
mergers of "be" , then we consider "abc" to be a neu- 
tralino candidate if /i = m(, c /m a (, c > A*min- The cut on 
z is used to avoid the soft splittings that dominate QCD 
branching and that are largely responsible for produc- 
ing the peak in the mass distribution of QCD jets. The 
cut on the subjet mass ratio fi ensures the presence of a 
three-body decay structure inside the jet. 

The full C/A- based analysis proceeds as follows: i) 
we use R — 0.7 and require at least four jets with 
p T > 500, 100, 100, 80 GeV, and \Ar n3 \,\A m4 \ < 1.5, 
ii) the hardest jet is taken to be a neutralino candidate 
if it passes the substructure cuts with z m ; n = 0.15 and 
Mmin = 0.25. Relative to the kx analysis, the harder 
jet cuts help keep the background under control despite 
requiring substructure in only one jet. For events that 
pass the cuts, we plot in Fig. [5] (right) the distribu- 
tion of m a b c , weighted with m a f, c /100 GeV. Expectations 
from QCD are that this distribution should be rather 
flat for m m i n < rrij < ptR\J 'Zmin, where m m i n is some 
small value governed by higher orders. This is indeed 



what we observe, and for a range of choices of R value 
and pt cut the signal is found to lie in this interval. 
The analysis reconstructs 10.3% of all neutralinos with 
Pt > 500 GeV (constituting 3.7% of all neutralinos pro- 
duced) in a 20 GeV mass window around the nominal 
mass; for QCD jets 0.052% of all jets with p T > 500 GeV 
are accepted in this mass window. 

For both analyses we estimate the significance of the 
signal from the number of signal and background events 
in a window of 20 GeV centered on the peak, a choice 
consistent with the 7% mass resolution seen in ATLAS 
detector simulations [l(J. This ignores the effect of 
the 'looking-elsewere' problem, but should demonstrate 
the potential of such a search. The results are shown 
in Fig. [3] (left) for various neutralino masses along the 
SPSla benchmark line. Even in the highest-mass case 
studied, the significance is well above the 5-tr discovery 
'threshold' with 1 fb" 1 of statistics. In fact, the order of 
magnitude loss in cross section over the parameter range 
is to a large extent compensated by the lower QCD back- 
ground at higher jet masses for the &t analysis, and in 
the C/A case, jets from higher-mass squark and gluino 
decays are more likely to pass our hard pt cut. 

With evidence for a resonance peak, the next step 
is to estimate the mass of the resonance. We fit the 
jet-mass distributions with a background plus Gaussian 
signal distribution. The fcy analysis uses an exponen- 
tial background in the interval [80, 200] GeV, while the 
C/A analysis uses a uniform background in the interval 
[75, 120] GeV. The results of this naive fit, which ignores 
the experimental jet mass resolution, are shown in Table|T] 
for the SPSla benchmark point. Improvements on the 
systematic errors inherent in this method are possible by 
calibrating the jet mass against the known masses of the 
W boson or top quark, for which a reasonably clean mea- 
surement should be possible in events in which one top 
quark decays leptonically. Improvements in mass mea- 
surement are also possible through filtering of jets [l2l ]. 

In Fig. [3] (right) we also demonstrate the potential of 
our method for reconstructing the squark mass. By se- 
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Analysis 




X 2 /ndf 


mg R 


X 2 /ndf 


fer, R=0.4 


99.0 ±0.5 


1.43 






C/A, R=0.5 


98.2 ±0.4 


0.85 


517.0 ±2.5 


1.27 


C/A, R=0.7 


98.4 ±0.6 


1.46 


526.6 ±3.0 


1.20 



TABLE I: Neutralino and squark mass fits for SPSla. The 



nominal masses are 



96.1 GeV and mg 



520 GeV. 



looting events from the C/A analysis in the signal band 
90 GeV < m a b c < 105 GeV and plotting both the mass 
of the neutralino combined with the third jet and com- 
bined with the fourth jet in the event, we arrive at the 
distribution in black, with a clear peak around 520 GeV. 
The interpretation of the peak is checked by plotting the 
sideband distributions, picking events from 75 GeV < 
m abc < 90 GeV (red) and 105 GeV < m abc < 120 GeV 
(green). These show no sign of a peak. By subtracting 
the averaged sidebands (dashed line), and fitting the re- 
maining peak with a Gaussian, we arrive at the squark 
mass estimates in Table HI 

The effects of pile- up, intrinsic resolution and granular- 
ity of the detector will all have additional impact on the 
discovery of a neutralino resonance and the measurement 



of its mass at the LHC, but initial studies with realis- 
tic detector simulations indicate that the efficiencies and 
resolutions assumed here are not unreasonable [lOi] . Fur- 
ther studies of the techniques presented here are ongoing 
within the ATLAS collaboration [Hj]. 

In conclusion, we see that using sophisticated jet clus- 
tering algorithms such as kr and C/A gives us the pos- 
sibility of discovering baryon-number violating decays of 
the type \i llli without the assumption of additonal 
features such as hard leptons, and even when using only 
the substructure of the hardest jet in the event. We have 
further found that the neutralino mass can be measured 
to a precision of a few GeV in these R-parity-violating 
scenarios, most likely limited by the experimental jet 
mass resolution, and that one can identify the squark 
resonance. Realizing the potential outlined in the above 
analyses is a challenge that merits experimental study. 
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